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Li and Haldane conjectured and numerically substantiated that the entanglement spectrum of 
the reduced density matrix of ground-states of time-reversal breaking topological phases (fractional 
quantum Hall states) contains information about the counting of their edge modes when the ground- 
state is cut in two spatially distinct regions and one of the regions is traced out. We analytically 
substantiate this conjecture for a series of FQH states defined as unique zero modes of pseudopo- 
tential Hamiltonians by finding a one to one map between the thermodynamic limit counting of two 
different entanglement spectra: the particle entanglement spectrum, whose counting of eigenvalues 
for each good quantum number is identical (up to accidental degeneracies) to the counting of bulk 
quasiholes, and the orbital entanglement spectrum (the Li-Haldane spectrum). As the particle en- 
tanglement spectrum is related to bulk quasihole physics and the orbital entanglement spectrum is 
related to edge physics, our map can be thought of as a mathematically sound microscopic descrip- 
tion of bulk-edge correspondence in entanglement spectra. By using a set of clustering operators 
which have their origin in conformal field theory (CFT) operator expansions, we show that the 
counting of the orbital entanglement spectrum eigenvalues in the thermodynamic limit must be 
identical to the counting of quasiholes in the bulk. The latter equals the counting of edge modes at 
a hard-wall boundary placed on the sample. Moreover, we show this to be true even for CFT states 
which are likely bulk gapless, such as the Gaffnian wavefunction. 

PACS numbers: 63.22.-m, 87.10.-e,63.20.Pw 



I. INTRODUCTION 

Determining the universality class of a real system ex- 
hibiting a topological phase is a difficult task in con- 
densed matter physics. Renormalization group meth- 
ods have been very successful in uncovering the universal 
physics in phases with local order parameters, but, due to 
their perturbative approach, cannot be readily general- 
ized to topological phases which do not exhibit symmetry 
breaking. The density matrix re nor malization group 1 -"^ 
and tensor matrix product states^ can probe topologi- 
cal order in one dimension, but have had limited success 
with higher dimensional systems so far. The prototype 
of two-dimensional topologically ordered phases are the 
experimentally accessible fractional quantum Hall (FQH) 
phases. A promising tool to extract topological informa- 
tion from the ground state wavefunction in these phases 
is the entanglement entropj^^. However, it depends on 
scaling arguments, is hard t o obta in to sufficient accuracy 
from numerical calculation d 12 ! 13 !, and does not uniquely 
determine the topological order in the state. 

In 2008, Li and Haldane 14 proposed a new tool to iden- 
tify topological order in non-abelian FQH states - the 
entanglement spectrum. They divided the single-particle 
orbitals in a Landau level on the sphere along the equa- 
tor and constructed the reduced density matrix of the 
ideal (model) and the realistic (Coulomb) FQH states in 
the upper half of the sphere (part ^4) by tracing out or- 
bitals in the lower half (part B). Having thus created 
a 'virtual' edge, they defined the orbital entanglement 
spectrum (OES) to be the plot of the negative logarithm 
of the eigenvalues of the reduced density matrix of A vs 
the z-angular momentum of A (Lf ) for a fixed number 



of particles in A. In particular, Li and Haldane consid- 
ered the part of the spectrum with the lowest-lying levels 
and the highest-weight eigenstates of the reduced den- 
sity matrix of A. They noticed that the number of levels 
in every OES of the model states, such as the Laughlin 
and the Moore-Read, was much smaller than the Hilbert 
space dimension and was identical to the counting of the 
conformal field theory (CFT) modes associated with the 
edge at small values of Lf. Although the number of levels 
in the OES of the Coulomb state saturated the Hilbert 
space dimension, a gap separated the levels higher in the 
spectrum from a CFT-like low-lying spectrum at small 
values of with the same counting as the model state. 
This was taken as evidence that the Coulomb state at 
v = 5/2 and the model Pfaffian state belonged to the 
same universality class. Based on extensive numerical 
evidence, they conjectured that — 1) In the thermody- 
namic limit, the counting of the OES of the model state 
is the counting of the modes of the conformal theory de- 
scribing its gapless edge excitations, 2) The 'entangle- 
ment gap' separating the low-lying, CFT-like levels from 
the generic ones higher in the Coulomb spectrum is finite 
in the thermodynamic limit. 

Many researchers have investigated properties of the 
entanglement spectra since. The authors of [15j discov- 
ered that the entanglement spectrum in the thin-annulus 
limit (the conformal limit) had for, several examples, a 
full gap at finite system sizes. The counting of the entire 
low-lying spectrum of the Coulomb state is the same as 
that of the corresponding model state in this limit. Mo- 
tivated by this result, we recently conjectured a counting 
principle for the finite-size counting of the OES of the 
Laughlin states^. Other cuts have also been studied. 
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FIG. 1. Left: Sketch of the partition in orbital space (part B in grey). The tracing procedure creates a virtual edge and the 
orbital entanglement spectrum (OES) probes the chiral edge mode(s) of part A. Right: OES of the u = 1/2 Laughlin state 
of N = 9 bosons, with orbital cut Ia = 8 and Na = 4. The minimal angular momentum L^ min defined in the text is the 
L^min — 20 sector in the plot. The entanglement level counting at = 20, 19, . . . , 16 is (1, 1, 2, 3, 5) which is the counting of 
modes of a U(l) boson in the thermodynamic limit. Finite size effects appear at = 15. 



Tracing out a fraction of the particles in the many-body 
ground state corresponds to a particle cut; the entangle- 
ment spectrum of the resulting reduced density matrix is 
the particle entanglement spectrum (PES) introduced in 
[17] . The level counting of the PES of a model state (de- 
scribed as a CFT correlator) is bounded from above by 
the number of bulk quasi-hole states of the model state; 
along with the OES, it is conjectured to contain all the 
topological numbers of the state. Entanglement spec- 
tra in other systems have also been explored; see, for 
instance, Refs. [T5rf55] . 

Analytic work in this emerging field is challenging be- 
cause of the strongly interacting nature of FQH states. 
The Li-Haldane conjecture (the correspondence between 
the counting of the number of modes in the thermody- 
namic limit OES and the counting of the edge-excitation 
spectrum) is easy to prove in non-interacting systems, 
such as the Integer Quantum Hall system and topolog- 
ical insulator a 36 " 38 !. In this article, we partially prove 
the first part of the Li-Haldane conjecture for cluster- 
ing model states: in the thermodynamic limit, we show 
that the counting of the low-lying levels of the OES is 
identical to the PES counting and is bounded above by 
the counting of the CFT associated with the edge. We 
prove this bound for the bosonic (fc, 2)-clustering states 
(the Read-Rezayi sequence) multiplied by any number 
of Jastrow factors, and for the Gaffnian. In principle, 
the bound should hold for all model states defined as 
the unique, highest-density zero modes of (k + l)-body 
pseudopotential Hamiltonians^l. The starting point of 
the proof is at the PES. We show that the counting of 
the PES is bounded above by the number of zero modes 
of the model pseudopotential Hamiltonian (number of 
bulk q-h excitations) of the particles in dr^. We then es- 
tablish a bulk-edge correspondence in the entanglement 
spectrum by showing that the particle and orbital en- 
tanglement spectra have the same counting. This cor- 
respondence holds for a part of every spectrum at finite 



size and for the entire spectrum in specific quantum num- 
ber sectors in the thermodynamic limit, as will be made 
precise later. Using the known bulk-edge correspondence 
in the energy spectrum of bulk gapped FQH states, we 
analytically prove that the counting of the natural OES 
is bounded above by the counting of the CFT describing 
the edge, and reason physically that the bound should be 
saturated. 

The paper is organized as follows: We define the or- 
bital and particle entanglement matrices and spectra in 
Sections [ITT] and [IV] In Sec. |IVC[ we present the upper 
bound to the number of levels in the particle entangle- 
ment spectrum and argue for its saturation. In Sec.[V] we 
formulate the clustering properties of the model state in 
the single-particle orbital basis and use them to relate the 
counting of the particle and orbital entanglement spectra 
of the Read-Rezayi sequence in Sec. |VI[ The proof for the 
upper bound of the Li-Haldane conjecture is presented at 
the end of the same section. In Sec. |VID[ we extend the 
proof to the other model states. Examples and the math- 
ematical formulation of the ideas in the proof are in the 
Appendices. 

II. NOTATION 

The results that we present in this article hold on any 
surface of genus (such as the disc or the sphere) pierced 
by Ntf, flux quanta; for simplicity, we choose the sphere 
geometry. The single-particle states of each Landau level 
are eigenstates of L Zl the z-component of angular mo- 
mentum and \L\ 2 , the square of the magnitude of the to- 
tal angular momentum vector^. In the Lowest Landau 
Level (LLL), the degenerate single-particle states belong 
to a multiplet of angular momentum L — and con- 

sequently, L z 6 [— . . . , iV^/2]. We are forced to 

adopt a dual notation in this article - the single-particle 
orbitals are indexed by their z-angular momentum L z in 
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FIG. 2. Left: Sketch of the partition in particle space. The particles of part B that are traced out are denoted in grey. 
Right: particle entanglement spectrum (PES) of the v = 1/2 Laughlin state of N = 9 bosons, with particle cut Na = 4. 
The minimum angular momentum L^ min defined in the text is L^ min = 20 in the plot. The entanglement level counting at 
Lz — 20, 19, . . . , 15 is (1, 1, 2, 3, 5, 6), identical to the counting of quasiholes in a Laughlin state of Na particles but at flux 
N$ = 16. 




FIG. 3. The configurations in the PES can be related to those of the OES using the clustering constraints. These constraints 
reveal the vanishing properties of the FQH state as particles are brought closer together. They relate the long-wavelength 
properties of the FQH state when two particles are far away from each other to the short-wavelength properties of the state 
when particles are close together, and hence can be used to 'drag' particles from the PES hilbert space into the more restrictive 
OES hilbert space. 



the figures, and by a shifted label — L z ) in the 

text. At the north (south) pole, L z = A^/2 {-N$/2). 

Fermionic/bosonic many-body wavef unctions of N 
particles and total angular momentum L f ° can be ex- 
pressed as linear combinations of Fock states in the occu- 
pation number basis of the single-particle orbitals. Each 
Fock state |A) can be labeled either by A or by n(A). 
A = [Ai, A2, . . . , Ajy] is an ordered partition of L| ot into 
N parts; each orbital with index Xj is occupied in the 
Fock state. By definition, Ai > Xj Hi < j. n(X) is 
the occupation number configuration. It is defined as 
n(X) = {rij(X),j = 0, ...N^}, where rij(X) is the occu- 
pation number of the single-particle orbital with angular 
momentum j. For instance, if = 6, the partition 
and occupation configuration of a Fock state |4, 2, 0) of 
3 particles is A = [4, 2, 0] and n(A) = {1010100} respec- 
tively. We will repeatedly run into a special kind of par- 
tition in this article — the (k, r)-admissible partition. A 
(k, r)-admissiblc partition labels a Fock state, whose oc- 
cupation configuration has no more than k particles in r 
consecutive orbitals. The partition in the example above 
is (1, 2)-admissible. 



Three useful relations between partitions are 'domi- 
nance', 'squeezing' and 'addition'. A set of partitions 
may always be partially ordered by dominance, indicated 
by the symbol '>'. A partition [i dominates another par- 
tition v (fi > v) iff J2 r i=0 Mi > £i= v i W e [°> • • • > N l 
Squeezing is a two-particle operation that connects n(p) 
to n(v). It modifies the orbitals occupied by any two 
particles in n(p) from mi and m-i to vn! x and m' 2 in n(v), 
such that mi +TO2 = m\ +ml 2 and mi < m\ < m' 2 < m 2 
if the particles are bosonic or mi < m'i < m! 2 < m 2 
if they are fermionic. Dominance and squeezing are 
identical concepts: a partition /z dominates a partition 
v iff v can be squeezed from /x by a series of squeez- 
ing operations. The 'sum' of two partitions /1 + v is 
defined as the partition with occupation configuration 
n{ji + v) = {nj{^) + n j {v),j = Q,...N ct> }. 

FQH wavefunctions in the LLL are translationally in- 
variant, symmetric, homogeneous polynomials of the co- 
ordinates of the N particles, (zi, z% . . . Zjy). We con- 
sider mainly the bosonic Read-Rezayi sequence at fill- 
ing v = fc/2 here (see Sec. VI D for other model 
states). These states are the unique, highest density 
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zero-mode wavefunctions of (k + l)-body pseudopoten- 
tial Hamiltonians^. Recent worlP-^ has shown the Read- 
Rezayi bosonic wavefunctions "0 to be Jack polynomials 
J" o indexed by a parameter a = — (fc+1) and the densest- 
possible (fc, 2)-admissible 'root' configuration^ 



n(X ) = {kOkOkO . . . kOk} . 



(1) 



For the (k, 2) states, the number of fluxes for the ground- 
state wavefunction is N<p = 2(N/k — 1). All the Jack 
polynomials at a = —(fc + 1) indexed by (fc, 2)-admissible 
root configurations are (fc, 2)-clustering polynomials, i.e. 



they vanish as ni>fe( z — z «) 2 when 



zi 



z k . They 



form a basis for all many-body (fc, 2)-clustering polyno- 
mials and can be decomposed into a linear combination 
of Fock states with configurations squeezed from the root 
partition. 



III. THE ORBITAL ENTANGLEMENT MATRIX 
(OEM) 

A. Definition 

Consider dividing the set of single-particle orbitals 
{0, 1, . . . , iV^} into two disjoint sets A = {0, 1, . . . Ia — 1} 
and B = {I a, ■ ■ ■ N<j,}. As the single-particle orbitals are 
polynomially localized in the 9 direction, this partition in 
the single-particle momentum space roughly corresponds 
to an azimuthally symmetric spatial cut. 

The number of orbitals in A (B) is I a (Ib), where Ib = 
N<f, + 1 — I a- Without loss of generality, let Ia < Ib 
(I a > Ib for A and B swapped). Any occupation number 
state | A) may be expressed as a tensor product Cgi \v) 
of states with partitions \x and v belonging to the Hilbert 
spaces of A and B respectively. Thus, the model state 
can be decomposed as: 



\ v i)i 



(2) 



i,3 



where the kets {l^i)} and {|^j}} form orthonormal bases 
that span the Hilbert spaces of A and B. The matrix 
Cf is the full orbital entanglement matrix (OEM). The 
(i, j)th matrix element of the full OEM is equal to the 
coefficient of + Uj) in |i/>), i.e. (Cf)y = b^ i+vr 

In this article, we will almost exclusively deal with en- 
tanglement matrices. Unless stated otherwise, the rows 
(columns) of these matrices, for both the OEM defined 
in Eq. ^ and for the PEM defined below, will be labeled 
by partitions /i.; (vj) corresponding to the occupation ba- 
sis states (\vj)) in A (B). The vector defined by the 
entries of a row/column in the entanglement matrix shall 
be referred to as row/column vector. 

The OEM is explicitly constructed for the 4-particle 
v = 1/2 Laughlin state with an orbital cut after 3 orbitals 
in Appendix | A 1| 



B. Properties 

Cf has a block-diagonal form; each block C is labelled 
by Na, the number of particles in A, and Lf, the to- 
tal z-angular momentum of the Na particles in A. Note 
that Lf = Y^i=i^i f° r the state where [li here are 
the components of the partition \l. Due to an unfor- 
tunate but necessary choice of notation, fj,j also index 
the partitions of the Hilbert space of part A. In that 
case, Hi is a partition by itself, and its components are 
Mii)Mi2i • • ■ j^iN A - The use of /ii as a partition or as a 
component of a partition fi will be self-evident in the text. 
To understand the origin of the block-diagonal structure 
of Cf, observe that is an eigenstate of the particle- 
number operator N and the total z-angular momentum 
operator L* ot . As both operators are sums of one-body 
operators, N = N a ®I+I®Nb and L\ ot = Lf (8)1+1® Lf. 
Thus, every |A) in Eq.Q is labelled by the quantum num- 
bers, N and L* ot , while every {\v 3 )) is labelled by Na 
(Nb) and (Lf). In the remainder of this article, the 
symbol C refers to the block of the full OEM Cf with 
labels Na and Lf. 

The reduced density matrices are obtained from Cf 
as pa = CfCf^ and pg = Cf^Cf. The block-diagonal 
structure of Cf carries over to the reduced density matri- 
ces and the rank of pA and ps in each block is equal to 
that of C. Neither pa nor pb uniquely determine all the 
coefficients of Cf clearly contains more information 
than either of the reduced density matrices. 

The singular value decomposition of C is given by: 



h3 



rank(C) 

E • 

i=l 



(3) 



The kets on the left-hand-side of Eq. (3j) are defined as 
in Eq. |Z7j) and \Vi) are the singular vectors in the 
Hilbert spaces of A and B restricted to a fixed particle 
number and z-angular momentum. They are linear com- 
binations of the occupation number basis vectors and 
\vj). The £i's are the 'energies' plotted as a function of 
Lf in the orbital entanglement spectrum (OES) intro- 
duced in [T3]. A typical OES is shown in Fig. [l] for the 9 
particle 1/2 Laughlin state with orbital cut I a = 8 in the 
sector Na — 4. Note the counting of the entanglement 
levels, {1, 1, 2, 3, 5 . . .}, starting from the right edge of the 
spectrum. 

The number of finite energies (rank(C)) at each Lf is 
independent of the geometry of the 2-d surface and the 
symmctrization factors arising due to multiple particles 
occupying the same orbital. Let C f d and C f s be the full 
OEMs in the disc and sphere geometry, or in any other 
two genus geometries. Modifying the geometry of the 
surface changes the normalization of the single-particle 
orbitals (the quantum mechanical normalization); thus 
every b\ in the expansion of \ip) in Eq. ^ in the disc basis 
is multiplied by a factor A/"(A) = n^i-^-^i)- when ex- 
panded in the single-particle orbital basis on the sphere. 
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Af(j) is a factor relating the normalization of orbital j on 
the disc to that on the sphere. The OEM's on the disc 
and the sphere are thus related as: 

i^>=x;(c f %iMf}®i^) 

= J2(c?%^(rfW(vfM) ® K), (4) 

where the superscripts d and s refer to the disc and 
sphere geometries, or to any other two genus ge- 
ometries. Cf S is obtained from Cf d by multiplying 
whole rows and columns by normalization factors; thus 
rank(Cf S ) = rank(Cf d ). An identical argument shows 
the rank of Cf to be independent of the symmetriza- 
tion factors that arise in the normalization of the many- 
body states constructed from normalized single-particle 
orbit als. We are therefore free to work in an unnormal- 
ized single-particle basis from this point. 

For a given cut I a in orbital space, the maximum num- 
ber of particles that can form a (k, 2)-clustering droplet 
in A is defined to be the natural number of particles 
N A , na t- Explicitly, N A , na t = k[(l A + 1)/2J where |xj 
is the integer part of x. The OES at N A nat is called 
the natural spectrum. In the thermodynamic limit, the 
counting of the OES is conjectured to be identical to the 
counting of the modes of the CFT describing the edge for 
values Ia,Na — > oo such that Ia/N^ — > const. (> 0) and 
N A /NA,nat -> 1. 

For future reference, L A min denotes the minimum z- 
angular momentum of the Na particles in A: 

Limin = V N A/k\ (2N A - k[N A /k\ - k). (5) 

We stress that L A min is the maximum value on the x- 
axis of the numerically generated entanglement spectra 
existing in the literature, due to the different indexing 
scheme in the text and the figures (see also the discussion 
in Section |n|. For instance, in Fig. [IJ L A min describes 
the sector of the OES at L A = 20. 



IV. THE PARTICLE ENTANGLEMENT 
MATRIX (PEM) 

A. Definition 

In the orbital cut that we just discussed, the Hilbert 
space of A at some (Na, L a ) was spanned by the possible 
occupation configurations \fi) of N A particles, such that 
Y]^i [ii = L A and fii < I a Vi. We now consider making a 
cut of a FQH state \ip) in particle space by dividing the N 
particles into groups A and B with Na and Nb = N—N A 
particles. Let N A < Nb- The Hilbert space of A(B) is 
now spanned by all possible occupation configurations of 
N a (Nb) particles in the full single-particle orbital basis 
of the state and contains the smaller Hilbert space of 



A(B) with the orbital restriction. Just as in the previous 
section, we can write the model wavefunction as: 

|^>=£&A|A)=£(Pf)dMi>®h->, (6) 

where the kets and are Fock states of Na 

and Nb particles in the full single-particle orbital space. 
The matrix P/ is the full particle entanglement matrix 
(PEM). As was the case for the OEM, the matrix ele- 
ments of the PEM are directly related to the weights of 
the model wavefunction by (P/)ij = b^ i+Vj . 



B. Properties 

For a given cut with N A particles in A, P/ is block- 
diagonal in the angular momentum of part A, The 
block of Pf at fixed (Lf,N A ) shall be denoted by P. 
The reduced density matrices of part A and B are given 
by pa = P/P/ an d Pb — P/P/ respectively. They are 
block-diagonal in and have the same rank as Pf in 
each block. In the same spirit as the discussion of the 
OEM, we define the singular value decomposition of the 
PEM by: 

X)(p)«k-) ® = £ e " £</2 i^> ® m, (7) 

i,j i 

where the singular vectors \Ui) and \Vi) are orthonormal 
vectors in the Hilbert spaces of A and B restricted to 
fixed angular momentum. The plot of the 'energies' vs 
is called the particle entanglement spectrum (PESjEI. 
In Fig. [2]we show the PES of the 9 particle 1/2 Laughlin 
state for the particle cut Na — 4. Note the same entan- 
glement level counting starting from the right edge of the 
spectrum, (1, 1, 2, 3, 5 . . .), in the OES in Fig. [l] and the 
PES in Fig. § 

In the spherical geometry, the PEM is labelled by 
an additional quantum number as compared to the 
OEM 17 — the total angular momentum of A, (L A ) 2 . 
Consequently, the eigenvalues of the block of the reduced 
density matrix with (L A ) 2 = £(£ + 1) have (21 + l)-fold 
degeneracy. This multiplet structure, apparent in the 
PES in Fig. [2] does not play any role in our discussions 
about the counting of the PES in this article. 

The OEM C at fixed particle number Na and angular 
momentum L A is a sub-matrix of the larger PEM, P, 
at the same angular momentum. L A min is the same for 
both cuts. 



C. Rank 

The property that defines the fc-clustered model state 
ip(z\, . . . , z/v) uniquely is that it is the lowest degree sym- 
metric polynomial that vanishes when (k+1) particles are 
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at the same position. Similar clustering conditions char- 
acterize every ground-state of a pseudopotential Hamil- 
tonian. This vanishing property must persist when we di- 
vide the particles into two groups and re-write the model 
state in Eq. ^ as: 



ip(zi, . . -,z N ) = 

^2 e ~ C * /2 ( z i> ■•■> z JVj \Ui) ® (zn a +i, z N \Vi), 



L A i 



(8) 



using Eq. Q at each Lf. If we choose [k + 1) particles 
in group A, say z%, . . . , Zk+i, to be at the same position 
z, then the state must vanish at every Lf. Further, as 
the singular vectors in B form an orthonormal basis: 



ip(z, z,z k+2 , ...,Z N ) = 



. Z, Zk+2, 



,z NA \Ui) = 0,Vi,Li. (9) 



A similar relation holds when A and B are interchanged. 
We conclude that the singular vectors, (zy, . . . , z^ A \Ui) 
and (zn a +i, ■ ■ ■ , z/v|Vj), must also be clustering polyno- 
mials that vanish when (fc + 1) particles are at the same 
position. A basis for clustering polynomials is the set 
of Jack polynomials, J~ , indexed by a = — (k + 1) and 
the (k, 2)-admissible partition jP^^. ip can therefore be 
expanded in the Jack basis as: 



i/}(zi, ...,z N ) 



■Zn a )J"Azn_ 



A + l: ■ ■ 



,z N ) 



(10) 



where fbi and Vj denote (k, 2)-admissible partitions of Na 
and Ng particles respectively. The matrix My is block- 
diagonal in angular momentum Lf; let M refer to the 
block of Mf at fixed value of Lf. The row and column 
dimensions of M are much smaller than those of P be- 
cause the (k, 2)-admissible partitions of Na and Nb form 
a small subset of the set of all partitions of with fixed Lf 
and Lf respectively. Nevertheless, as Eq. (|8| and (10) 
are equal, M and P must have the same rank. Thus, the 
minimum of the row and column dimension of M bounds 
the rank of the PEM from above. 

Let us reformulate what we have just shown in a more 
familiar language and argue for the saturation of the 
bound. The smaller dimension of M is the row dimen- 
sion because Na < Nb- It is equal to the number of dis- 
tinct bulk quasi-hole excitations of the (fc, 2)-clustering 
model state of Na particles at angular momentum Lf 
on a sphere pierced by the number of fluxes of the origi- 
nal state, = 2/k(N - k). For Lf < Lf mm + [N A /k\ , 
this is the number of distinct bulk quasi-hole excitations 
in the thermodynamic limit, i.e. no finite-size effects oc- 
cur. We have thus bounded the level counting of the PES 
by the number of bulk quasihole excitations at each Lf. 



Without further symmetry-induced constraints on the re- 
duced density matrices (we have already used all the sym- 
metries available in the state), we expect this bound to 
be saturated. In the thermodynamic limit (N A ,N — > oo 
such that Na/N > 0), we therefore argue that the level 
counting of the entire PES is identical to the number of 
the bulk quasi-hole excitations. This bound saturation 
can be proved exactly for the Laughlin states^. 

It is beneficial to identify a set of rows and columns 
in P with the same rank as the full matrix. Consider 
the rows and columns labelled by the (fc, 2)-admissible 
partitions. This sub-matrix of P is denoted by P and 
has the same dimensions as M. In Appendix |B 1| we 
show that P and M have the same rank. P will play 
a prominent role in the proof establishing the bulk-edge 
correspondence in the entanglement spectra. 

Configurations of Nb bosons with Lf 
< > 



O ^ 



bjj 




O 

O 



FIG. 4. A cartoon of the various sub-matrices in P 



V. CLUSTERING CONSTRAINTS 

In this section, we introduce the (&+l)-body clustering 
constraints that relate the rank of the PEM and the OEM 
of the clustering model states and establish the bulk-edge 
correspondence in the entanglement spectra. The Read- 
Rezayi model wavefunctions tp(k,2) i z i, ^2 ■ • • zn) are sin- 
gle Jack polynomials labelled by a root partition Ao 
(Eq. ([IJ), and a parameter a = —{k + 1). They sat- 
isfy (k, 2)-clustering — they are non-zero when a cluster 
of k particles is at the same point in space z = Z\ = 
Z2 = . . . Zfc , but vanish as the second power of the dis- 
tance between the (k + l)st particle and the cluster as 
z k+i z. The clustering property imposes a rich struc- 
ture on ip(k,2)( z i, Z2, ■ ■ ■ , zn). All the partitions A that 
arise in the expansion of l^) in the many-body occupa- 
tion basis = ^2\b\\X)) are dominated by Ao- Fur- 
thermore, all the coefficients b\ are known up to a mul- 
tiplicative constant. In the Jacks, this constant is chosen 
so that b\ a = 1. In other words, the clustering property 
and the requirement to be the densest possible wavefunc- 
tion determine ip(k,2) ( z i, Z2, ■ ■ ■ , z^) uniquely up to an 
overall normalization constant. Here, we formulate the 
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conditions imposed by clustering on tj}(k,2}{z\, . . . , zn) as 
linear, homogeneous equations on the coefficients b\. 



A. Derivation 

Let us introduce a 'deletion' operator for orbital i 
such that: 



di\X) = 







|a\{*» , 8 a 



(ii) 



A\{z} is the partition with a single occurrence of the or- 
bital i removed from it. The 'deletion' operators com- 
mute with each other. In Appendix [C] we derive the 
relation between these operators and the annihilation op- 
erators in the normalized single-particle basis. 

We now separate the coordinates of fc + 1 particles from 
the rest and rewrite ip(k,2)( z ii z 2. 




z k+2, 



, Zn) as: 



fc+i 



(12) 



and form a cluster by bringing the fc particles with coordi- 
nates z±, . . . , Zfc to the same position z. When z^ + i = z, 
the LHS vanishes and Eq.(12) becomes: 



E • 

h,...,l h+1 =0 



Y- k + 1 I ■ i 

^ = 1 3 (Zk+2,- 



fc+1 

,Zn\ JJ d L . 

i=l 



(13) 



The right-hand-side is a polynomial in an arbitrary com- 
plex number z, and has to vanish for every power f3 — 



Ylj—i h 01 z to satisfy the above equation. Thus, the 
constraints on \^) are: 

fc 

E d P-T.U h II d h ] M = D pWi = °- (I 4 ) 



ft is the z-angular momentum of (fc + l)-particles; it 
ranges from to N<f,(k + 1). The equation above re- 
quires any clustering wavefunction to be simultane- 
ously annihilated by the destruction operators {Dj , i = 
Q...Nt(k + l)}. 



The set of linear, homogeneous equations in Eq.(14| 
are linearly dependent. The dimension of the null-space 
of the set is exactly one for the densest possible wavefunc- 
tion, i.e. the vector of coefficients {b\} is uniquely deter- 
mined up to an overall multiplicative factor. Since the 
solution to Eq. ( 14 ) causes ip to vanish when any cluster 



of size greater than fc is formed in real space, we conclude 
that the set in Eq. ( 14 ) includes all constraints imposed 
on ip{zi, ■ ■ ■ , zn) due to clustering. 

Equivalently, we are describing model FQH wavefunc- 
tions that are the unique, highest density zero-modes of 
the Haldane pseudopotentials or their generalization to 
the fc + 1-body interaction In fact, the destruc- 

tion operators above are the fundamental clustering op- 
erators from which the Haldane pseudopotentials can be 
obtained as the translationally-invariant supersymmetric 
form: 



(15) 



f(j3) can be derived at each fc; in Appendix [C] we work 
through the k = 1 case. 



VI. RELATING THE OES AND PES 
COUNTING 

We now have all the ingredients necessary to relate the 
counting of the PES to that of the OES for a given num- 
ber of particles in A and cut in orbital space. In Sec. |IV C| 
we constructed the sub- matrix P of P, labelled by (fc, 2)- 
admissible configurations with the same rank as P. We 
now use the clustering relations derived in the previous 
section to express the row/column vectors of P in terms 
of those of the OEM C. When possible in finite-size, 
i.e. up to a certain value of Lf , this procedure proves 
that the PES and the OES have the same counting. In 
the thermodynamic limit, this procedure establishes the 
equality of the level counting of the entire PES to OES 
when NA/NA,nat L thus proving a significant part of 
the Li-Haldane conjecture. 

The argument below applies equally well to row and 
column vectors. To keep the discussion concise, we for- 
mulate it using row vectors alone. 



B. Properties 



Every value of /? in Eq. (14) yields, in general, a large 



number of linear relations between the coefficients of \ip). 
Let Sp be the set of all partitions of N particles such that 
the sum of the z-angular momentum of (fc + 1) particles 
is (3. For every occupation configuration of N — (fc + 1) 
particles, Eq. ( 14 1 relates the coefficients of partitions A € 
Sp in the expansion of \ip). Examples of such relations 
are given in Appendix [P] 



A. Systemizing the constraints 

The biggest challenge in relating the row vectors of 
P to those in the OEM lies in identifying a set of lin- 
early independent equations in the entire set of clustering 
constraints. To this end, we introduce a few quantities 
characterizing a partition /i. n m (/i) below refers to the 
occupation number of the with orbital in partition u. The 
orbital cut is after I a orbitals. 



The unit cell — We divide the single-particle orbital 
space such that the jth unit cell contains the orbitals of z- 
angular momentum 2j and 2j + l, and j € [0, . . . , N^/2). 
As the total number of single-particle orbitals is odd for 
the bosonic (k, 2)-clustering states, the orbital with an- 
gular momentum is its own unit cell with index N^/2. 
Every orbital belongs to exactly one unit cell. 

The intact unit cell — The jth unit cell of a partition 
/x is said to be intact if the occupation numbers of the 
orbitals with angular momentum 0, . . . , 2j + 1 arc iden- 
tical to those in the root configuration Eq. ([!]), i.e. if 
n-i (/i) = nj(Ao) for i = 0, 2j + 1. Clearly, the jth unit 
cell can only be intact if all unit cells 0, . . . , j — 1 are 
intact. 

The number of intact unit cells in part A — The num- 
ber of intact unit cells in part A, A M , is the number of 
intact unit cells to the left of the orbital cut in n(/z). 

Distance from the cut — If we were to number the or- 
bitals to the right of the cut as 1, 2, . . ., then the distance 
from the cut is defined as the sum of the indices of the 
occupied orbitals to the right of the orbital cut in n(n). 
The distance from the cut, K^, is given by: 



Ky. = k + 

12 3 4 


3, 1 a = 1C 
5 6 


kO 


kO 


kO 


0k 


33 


kO 


1 


■* ► 

A M = 3 

12 3 4 


5 6 


k0 


kO 


kO 


0k 


313 


kO 


1 



A M = 3 



K„ = 2k + 7 
I A = 9 



FIG. 5. The occupation configuration of a generic partition fj, 
with the unit cells, the number of intact unit cells A M and the 
distance from cuts after I a = 9 (top) and Ia = 10 (bottom) 
shown. N,i, = 12 here. 



The method 



Consider the blocks P and C at Na and Lf = Lf min + 

£. We can relate all row vectors of P to row vectors in 
the OEM if: 



A Ao >K fio , 



(17) 



K„ 



m—lA 



n m (p)(m ~l A + 1). 



(16) 



where n(fio) is a (fc, 2)-admissible occupation configura- 
tion and labels a row of P: 



n(/I )={ fcO.^.fcO xO^OlO.-.O}. 

2L(i\U-l)/fcJ i-\ 



(18) 



K(fi) — for a partition /i labeling a row of the OEM; for 
a general partition, it represents the distance in orbital 
units that all the particles to the right of the cut need to 
traverse to cross the cut. In Fig.[5j we pick as an example 
a generic partition \i and identify the number of intact 
unit cells in A, A M , and the distance from the cut, K^, 
for two different orbital cuts. 



The value of x is fixed by the total particle number being 
N A (x = (N A - 1) - k[(N A - l)/k\). /i dominates all 
the other partitions that label rows of the PEM. It has 
the maximum (total) number of intact unit cells possible, 
l(N A — l)/fcj . The constraint, A^ > Kjx , sets an upper 
bound on the values of I for which we can express all row 
vectors in P in terms of those in C: 



A Ao - fcA? o - (2A Ao + l)(a: + 1) for l A even, l A < 2[(N A - l)/k\ 
<{ A Ao - fcA^ (A Ao - 1) - 2{x + l)A Ao for l A odd, l A < 2[(N A - l)/k\ (19) 
I A - A Ao + 1 for I A > 2 [(N A -l)/k\, 



where we abbreviated the difference of the total number 
of intact unit cells and those only in A by A^ = [(-^Va — 
l)/fcj — A^ . For such values of £, the method of proof 
can be broken into two steps — 

I If A^ > if/j , then A^ > for all (fc, 2)-admissible 
partitions fi < fiQ. 

II If A^ > for a partition /z, then the row vector 



labelled by /j in P can be expressed as a linear com- 
bination of row vectors in the OEM C alone. 

We prove these statements rigorously in the Appen- 
dices [E] and [F] The first step shows that the A^ > Kjx 
for all partitions fi labeling rows of P; the second assures 
that all these rows can be written as linear combinations 
of rows in the OEM alone. 
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An identical argument can be repeated for the 
columns; let vq be the (k, 2)-admissible partition that 
dominates all partitions of L? into Ng parts. For values 



of I such that A^ > Kp_ and Aj 



> Kc , the OEM 



and PEM have the same counting in finite-size. I is usu- 
ally of order N^/k. We can therefore analytically prove 
the bulk-edge correspondence in the entanglement spec- 
tra for a few values of Lf near L^ min (the right edge of 
the numerically generated spectra). Note that at these 
values of £, it was proved in [17] that the counting of the 
PES is bounded from above by the number of bulk quasi- 
hole excitations in the thermodynamic limit. We argued 



in Sec. IV C that due to the absence of extra-symmetries, 
this bound should be saturated. 

The heart of the proof lies in the use of the (k + 1)- 
clustering condition at the z-angular momentum of the k 
particles in the right-most intact unit cell in part A and 
one particle occupying an orbital to the right of the cut. 
This relates a single row vector belonging to the PEM 



with A M and to row vectors with A M / 



Ap — 1 and 



Ku' < Ku ~ 1- The clustering conditions thus allow us 
to replace a row vector whose partition has distance 
with a linear combination of row vectors whose partitions 
have distances reduced by at least one at the cost of using 
a single intact unit cell. If A p > for a partition //, 
then iterating this procedure provides a linear relation 
between the row vector labelled by partition /i and row 
vectors with distance zero, i.e. row vectors of the OEM 
C. 

In the thermodynamic limit, Na and Ia scale with TV. 
Consequently, the number of intact unit cells in A (73) 
in /Sq (£o) denoted by A^ [Aj> ) scales with TV. For the 
upper bound on £ in Eq. ( 19 ) to scale with TV, A^ and 

TV: 



Ap have to grow slower than 



A^Tq , Ayg 



A M ~ , A;7 



TV 

N 9 , g < 1/2 . 



As: 



\N A -N A ,nat\ 



\N B - N B . 



\Na — NA,nat\ must grow slower than vTV. Thus, in 



the thermodynamic limit, condition (19) can only be 
fulfilled if NA/NA,nat — ► 1- In OES sectors such that 
NA/NA,nat — > 1, the level counting of the OES is then 
identical to the level counting of the PES at Na — Na , n at 
for all angular momenta L^, thus proving the bulk-edge 
correspondence in the entanglement spectra. 



C. Illustrative examples 



1. Atk = l: 

Consider the v = 1/2 Laughlin state of TV = 7 bosons 
with TV = 12 and L* ot = 42. Let l A = 6 and the number 
of particles in A be the natural number Na = NA,nat — 3. 
We consider the entanglement level counting of the OES 



and the PES at Li 



L z , min + £ = L z<min + 3. We first 
verify that the conditions, A^ > K^ and Ap > K„ , 
are satisfied. The occupation configurations of /Jo and isq 
are: 



n(/T ) = {101000|0100000} 
n (z7 ) = {000100 | 0010101} 



Ku 



: 2,A^ 
3, A* = 



= 2 
3 



The cut in orbital space is indicated in the occupation 
configurations by the '|' symbol. Hence, the method dis- 
cussed in the previous section should prove the equality 
of the ranks of the OEM and the PEM at this Lf. 

The occupation configurations of the (1, 2)-admissible 
partitions labeling the rows of P are: 



n(iXo) = { 101000 | 010... 0} 
= {100100| 100... 0} 
n(jl 2 ) = {010101 1 000... 0} 



= 0, A^ = . 



(20) 



/12 labels a row that already belongs to the OEM C. We 
now relate the row labeled by the partition fix to rows of 
the OEM. In n(/Ji), only the 0th unit cell is intact and the 
particle to the right of the cut occupies the orbital with 
index 6. Let us pick the 2-body clustering constraint at 
(8 = 6 (the sum of the z-angular moment A of the particle 
in the intact unit cell and the particle to the right of the 
cut) in Eq. fify: 



{2(d d 6 + dxd 5 + d 2 di) + d 3 d 3 )\ip) = (21) 

For every occupation number configuration of (TV — 2) 
bosons with angular momentum (L^ — j3), Eq (21 ) gives 



one linear relation. The appropriate occupation number 
configuration for our purpose is n([3] + Uj), as: 



eMe {\fa + Vj)) = | [3] + Vj) . 



(22) 



into Nb = 4 parts label the 



The partitions Vj of Lf 
columns of the PEM P. 
indexed by fa to row vectors indexed by following parti- 
tions: 



Eq. (21) then relates the row 



= {010101 1 000... 0} 
n(jj, 2 ) = {001110 | 000... 0} 
n(jj, 3 ) = {000300 1 000 ... 0} 



=0,A Ml =0 



K„ 



A" 



0,A A 
0,A L 



(23) 



The proof of the full method is presented in the ap- 
pendices; here we illustrate the more formal ideas with 
examples of the general method at work for the k = 1 , 2 
wavef unctions. 



At every column index j, the explicit relation from 
Eq. plh is: 



2(P 



2jJ 



Paj = 0, 



(24) 
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where Py is the coefficient in P of the row labeled by /Xj 
and column labeled by Vj. We have thus related a row 
indexed by a partition fli with = 1 and = 1 to 
rows indexed by partitions /ii,/i2,/i3 with distance from 
the cut reduced by 1 and number of intact unit cells in A 
reduced by 1. These partitions label rows in the OEM in 
this example. A similar procedure, using the additional 
clustering constraint at f3 = 9, involving the particles in 
the orbitals of angular momenta 2 and 7, can be used to 
relate the row of P indexed by the partition /2o to rows 
in the OEM. 



2. Atk = 2: 

Let us now consider the Moore-Read state with N = 
18, N<f, = 16, and = 144 and perform an orbital cut 
after Ia = 7 orbitals. Here, we are interested in the sub- 
block of P and C with N A = 8, Lf = Lf mm + i = 

Lf min + 3. The occupation number configurations of /2q 



(2) 

where contains all terms involving angular momen- 
tum orbitals and/or 1. 

The clustering conditions in Eq. ( 26 ) yield a linear re- 
lation between certain coefficients in for each occu- 
pation number configuration of the remaining N — 3 par- 
ticles. We choose the configurations n([7, 5, 4, 0, 0] +Vj), 
as: 

|[7,5,4,0,0] + i/ J ) = d 2 (fadr(lA*3 + ^)), (28) 

where the \ vj) label the column vectors of the PEM. Note 

that d 2 d 2 d^ is the only term in that contains the an- 
gular momentum 7 orbital; all other terms have highest 
angular momentum less or equal 6, and thus smaller dis- 
tance to the cut. Equivalently we can note that as D\ x 
annihilates any configuration with an occupied orbital 
of z-angular momentum greater than 7, the first term 
in Eq. (26) relates the row labeled by ^3 only to rows 



labelled by partitions that are dominated by /I3: 



and vq are: 






= {2011111 


1000. 


.0} 




= 1,A M1 


= 1 


n(jl ) = {2020201 1 00100000} 


= 3, Afi = 3 


n(p 2 ) 


= {2010220 


1000. 


.0} 




= l.A w 


= 1 


n(u ) = {0000010 | 01020202} K u ~ 


= 2,A v - =3, 


nips) 


= {2002120 


1000. 


.0} 




= 1>A M3 


= 1 


where we indicate the orbital cut by the ' 


' symbol. Thus, 


n(fi 4 ) 


= {2001310 


1000. 


.0} 




= 1,A M4 


= 1 



and Aa 



> Kf) a , and we can relate all rows 



and columns of the PEM to ones in the OEM. 

The occupation configurations of the (2, 2)-admissible 



partitions labeling the rows 



of P are given by: 



n(fl ) 


= {2020201 


0010. 


.0} 




= 3, A M ~ 


= 3 




= {2020200 


1100. 


.0} 




= 3,A^ 


= 3 


"(Ma) 


= {2020111 


0100. 


.0} 




= 2, A„ 2 


= 2 




= {2020110 


2000. 


.0} 




= 2,A A 


= 2 


n{fii) 


= {2011111 


1000. 


.0} 




= 1,A„ 4 


= 1. 
(25) 



(29) 

All the partitions above have one less intact unit cell, and 
smaller distance K^. = Kp, 3 — 1 from the cut as compared 
to /Et 3 . 

The second operator in the clustering condition 



Eq. ( 26 1 acts on states with occupation number configu- 



rations such as: 



The trailing 0's in every occupation configuration indi- 
cate that the orbitals with L z = 10, . . . , 16 are unoccu- 
pied in the partitions labeling the rows of P. > Ka 
is satisfied for all i = 0, . . . , 4, as required in step [T] in 
Sec. rVTBl 

We illustrate the use of the 3— body clustering con- 
ditions by relating the row labelled by the partition jl 3 
to rows labelled by partitions (Xj with distance K^. = 1 
from the cut. The first unit cell is the rightmost intact 
unit cell in A in n{pi 3 ). Consider the 3-body clustering 
condition at j3 equal to the z-angular momentum of the 
2 particles in the rightmost intact unit cell and a particle 
to the right of the cut, i.e. at /3 = 11 = 2 x 2 + 7. It is 



beneficial to divide the clustering condition ( 14 1 into two 
terms: 



{4000110 1 100010 ... 0} 
{3100110 1 100100... 0} 
{3010110 1 101000... 0} 
{3001110 1 110000... 0}. 



All the above configurations have distance from the cut 
larger than K^ 3 = 2, and more than 2 particles in an- 
gular momentum orbitals and 1. Hence, they are not 
dominated by the root partition Ao, and have zero weight 
in the model wavefunction (the corresponding row in the 
PEM is identically 0). 

Thus, the clustering condition at /3 = 11 for the config- 
uration of the remaining particles being n([7, 5, 4, 0, 0] + 
i>j), yields a linear relation between the row labeled by 
/} 3 and the rows labeled by the partitions fli, . . . , /j, 4 : 



K ) +<)i^> 







(26) 



3 3 



2P 



2P 



2j 



3j 



*3 



0. 



(30) 



D$ = d 2 d 2 df + 2d 2 d 3 de + 2d 2 d 4 d$ + d^d^d^ + d 3 d 4 d 4 

-,(2) 



-Dii = d d dii + 2d d 1 d 10 + 2d d 2 d 9 



(27) 



where Py is the coefficient in P in the row labeled by /i.; 
and column labeled by Vj . The rows labeled by fj,\ , . . . , pb^ 
can in turn be related to rows in the OEM by using the 
clustering conditions at /? = 7. 
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D. Beyond (k, 2)-clustering states 

Until now, we have restricted our discussions to the 
bosonic (k, 2)-clustcring states ip(k,2){ z i, ■ ■ ■ , z n)- In this 
section, we generalize our results to other states with the 
property of clustering - the states obtained by multiply- 
ing (k, 2)-clustering states with M Jastrow factors and 
the (2, 3)-clustering Gaffnian state. We believe that our 
results hold for all highest-density states uniquely de- 
fined by clustering, like, for instance, the Haffnian state. 
For th e non -unitary states, which are supposedly bulk 
gaplesPM] the map relates the counting of the OES to 
the number of bulk quasihole states (which is equal to the 
counting of the PES); however, in this case the number 
of bulk quasiholes is not equal to the number of the edge 
modes, as the edge-bulk correspondence in the energy 
spectrum does not hold for non- unitary states. 

For the (k, 2)-clustering states, we identified a sub- 
matrix of the PEM, P, with the same rank as the PEM 
and whose smaller dimension was the number of distinct 
bulk quasi-hole excitations. We then argued, based on 
the lack of other symmetries in P, that its rank was equal 
to the smaller dimension, and that the PES counted the 
number of bulk quasi-hole excitations at each angular 
momentum. To generalize this argument to other clus- 
tering states, we need to first identify the special sub- 
matrix P. We can then establish the bulk-edge corre- 
spondence in their entanglement spectra by slightly mod- 
ifying the method used in Section |VI B| Extending the 
ideas in Sec. [VI] is quite straightforward - we re-define 
the notion of unit cell and the intact unit cell, and iden- 
tify N c , the number of linearly independent clustering 
conditions that involve the k particles of an intact unit 
cell and one particle to the right of the cut, for a fixed 
occupation configuration of the remaining N — (k + 1) 
particles. N c — 1 for the (k, 2) case. Using the iV c inde- 
pendent linear equations, we can relate a row labelled by 
a partition fi with A M intact unit cells and distance to 
the cut to rows labelled by partitions //, such that 
A pi = An — 1 and < — N c . Thus, in the notation 
of Sec.fVll when > K fio /N c and A v - > K„~JN C , the 



consecutive orbitals. The unit cell has (Mk + 2) or- 
bitals and the occupation configuration of the intact unit 



OES and the PES have the same counting. In the ther- 
modynamic limit, the arguments in the last paragraph in 
Sec. |VI| show that the Li-Haldane conjecture is true for 



these states as well. 



1. The (k, 2) -clustering state multiplied by Jastrow factors 



Let us consider the model wavefunction: 
i/)(zi, ...z N )= i>( k ,2)(zi, ■■■zn) Y[( z i ~ z j) M ■ 



(31) 



where ip(k ,2)( z i, ■ ■ ■ z n) is the (k, 2)-clustering state. In 
Appendix 



B2 



we show that P is labelled by row and col- 
umn occupation configurations that obey the generalized 
Pauli principle: no more than 1 particle in M consec- 
utive orbitals and no more than k particles in Mk + 2 



cellis^O)^ 1 !^) 



\Af-l 



1(0) M " 1 00} with^O)*'- 1 re- 



peated k times (we could succintly write the whole pat- 
tern {(l(0) Jl/_1 ) fe 00}). The exponent is the number of 
times the pattern in the parenthesis is repeated. In Ap- 
pendix [G^J we show that N c = 1 for M = 1. More 
generally, N c = [M/2\ +1 for k = 1 Laughlin states, and 
N c = 2[M/2J + 1 for states with k > 1. 



2. The Gaffnian state 

The Gaffnian state is a (2, 3)-clustering state and is a 
single Jack polynomial: 



ip(zi, . . . z N ) = Jx (zi, • ■ ■ z N ) 



(32) 



where a = -3/2 and n(A ) = {200200 . . . 2002}. It isde- 
scribed by a non-unitary CFT, the W 2 (3, 5) modePHU. 
It has been suggested that the fermionic Gaffnian state 
is the critical state between a strong-pairing phase and 
a Read-Rezayi phased Despite the Gaffnian being a 
gapless state, we can determine the counting of the PES 
and establish the correspondence in counting between the 
orbital and particle entanglement spectrum. The dis- 
cussion in Sec. |IV C| and Appendix |B 1| applies to any 
Jack polynomial with (fc, r) clustering that is a unique 
zero mode of a pseudopotential Hamiltonian (besides the 
(k, 2) Jacks, only one other Jack (2, 3) - the Gaffnian - 
satisfies this constraint). P is therefore the sub-matrix 
of the PEM labelled by (2, 3)-admissible row and col- 
umn occupation number configurations for the Gaffnian 
state. The unit cell has 3 orbitals and the occupation 
configuration of the intact unit cell is {200}. We derive 
the clustering conditions in Appendix |G l| and show that 
N r = 2 for the Gaffnian. 



3. The Haffnian state 

The Haffnian stateP is a (2, 4)-clustering state, but is 
not a single Jack polynomial. We cannot rigorously iden- 
tify P for the Haffnian state, although we expect, based 
on our understanding of the other model states, that P 
only contains the rows and columns labeled by partitions 
obeying the generalized Pauli principle discussed in Ref. 
|51) . We have verified this numerically. The occupation 
configuration of the intact unit cell is {2000}. The clus- 
tering conditions arc derived along the same lines as for 
the Gaffnian in Appendix |G 1[ giving N c = 3 for the 
Haffnian. Whenever > Kp; /3 and A„ > we 
numerically observe that the ranks of the PEM and the 
OEM are equal. 
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VII. CONCLUSIONS 

In this paper we have provided a proof that the Li and 
Haldane natural entanglement spectrum in the thermo- 
dynamic limit is bounded from above by the number of 
modes of the CFT describing the edge physics. Barring 
the presence of extra accidental symmetries in the sys- 
tem, the bound should be saturated. We showed that a 
part of the levels of two different entanglement spectra, 
the OES and the PES, probing different physics - that 
of the edge and bulk respectively - can be mapped into 
one another by a series of newly introduced clustering 
operators. This map works, for certain quantum num- 
ber sectors, both in finite size and in the thermodynamic 
limit and provides a mathematically sound foundation for 
a edge-bulk correspondence in the entanglement spectra. 
Our map works for both unitary and non-unitary states 
that are defined as unique highest density zero-modes of 
Haldane pseudopotential Hamiltonians. 
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Appendix A: A simple example 



the bulk-boundary correspondence in the following sub- 
sections. 



1. The orbital cut 

Let us cut the single-particle orbital space after I a = 3 
orbitals. Consider the blocks of the OEM at the natu- 
ral number of particles in A, Na = NA,nat — 2. From 
the above decomposition, the minimum possible angular 
momentum, Eq. |5]) for 2 particles in A is Lf min = 2. 
At this Na and Lf, the Hilbert spaces of A and B are 
spanned by = |2,0), |/x 2 ) = |1, 1) and \v x ) = |6,4), 
W2) = |5,5) respectively. The block C at Na = 2 and 
= 2 is then given by: 



|6,4) |5,5) 
|2,0) ( 1 -2 
11,1) \ -2 4 



(A2) 



where we have indicated the states labeling the rows and 
columns. Cy = b^+v. ( '+' as defined in Sectionpl|) and 
rank(C)=l. ' 

The block C with N A = 2,Lf = Lf mm + 1 = 3, of 
rank 1, is: 



|6,3) |5,4) 
|2,1) (2 -2 ) 



(A3) 



The block C at Na = 2, = Lf mm + 2 = 4, also of 
rank 1, is given by: 



|5,3) |4,4) 
|2,2) (4 4 



(A4) 



Fig. shows the numerically generated OES for the 
4 particle Laughlin state in the sphere geometry at 1/2 
filling with Na = 2 and I a = 3. The counting of the 
entanglement levels in the spectrum equals the ranks of 
C at each Lj-. 



Let us consider the bosonic Laughlin wavefunction of 
N = 4 particles at filling v = 1/2. The number of flux 
quanta, N$ is 6 and L' ot = 12. The wavefunction \ijj) can 
be expanded in the unnormalized basis as: 

|V> = £&a|a) 

A 

= 16,4,2,0) -2|6,4, 1,1) -2|5,5,2,0)+4|5,5,1,1) 
+ 2|6, 3, 2, 1) - 2|5, 4, 2, 1) + 4|5, 3, 2, 2) + 4|4, 4, 2, 2) 

- 2|6, 3, 3, 0) + 2|5, 4, 3, 0) - 6|4, 4, 4, 0) - 4|5, 3, 3, 1) 

- 6|6, 2, 2, 2) + 4|4, 4, 3, 1) - 6|4, 3, 3, 2) + 24|3, 3, 3, 3). 

(Al) 

We construct several orbital and particle entanglement 
matrices, and use the clustering constraints to prove 



2. The particle cut 

Let us construct the entanglement matrices for the par- 
ticle cut with Na = 2. At the smallest possible angular 
momentum = Lf min = 2, the PEM and OEM are 
identical: 

|6,4) |5,5) 

!tij ( -2 4 ) • (as) 

The Hilbert space of A at = L^ mm + 1 = 3 is 
spanned by the occupation number states |3, 0) and |2, 1). 
1 3, 0) was not a member of the Hilbert space of A for the 
orbital cut after I a = 3 orbitals (discussed in the previous 
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section), because the orbital with index 3 belonged to B. 
The PEM at Lf = 3 is given by: 



|3,0) 
12,1) 



|6,3) 
-2 
2 



M) 

2 

-2 



(A6) 



We see that the OEM (A3 1 for N A = 2, Lf = 2 is indeed 



a sub-matrix of the PEM, as discussed in Section IV A 

As the last example, consider Lf = Lf min + 2 = <±. 
The row and the column dimension of the PEM is larger 
than that of the OEM in (|A4|: 



(A7) 





J6,2) 


|5,3) 


|4,4) 


4,0) / 




2 


-6 


3,1) 




-4 


4 


2,2) \ 


v " 6 


4 


4 



By inspection, we see that the OEM (A4) is the sub- 
matrix consisting only of the first row and the first two 
columns. The rank of the PEM at Lf = 4 is equal to 
two and greater than that of the corresponding OEM. 

Fig. |6jb) shows the numerically generated PES for the 
4 particle 1/2 Laughlin state for Na = 2. The counting 
of the spectrum agrees with the ranks calculated above. 




(b) 



FIG. 6. (a) Orbital entanglement spectrum of the v — 1/2 
Laughlin state with N = 4, N A = 2, and orbital cut after 
Ia — 3 orbitals. The entanglement level counting is equal to 
the rank of the OEM at each angular momentum, (b) Particle 
entanglement spectrum of the v = 1/2 Laughlin state with 
particle cut Na = 2. The entanglement level counting at all 
angular momenta Lf is equal to the rank of the PEM. L 
defined in the text is Lf = 4 in the plots. 



3. Relating the OES and PES counting 



A 
z.min 



At Lf = 2, C (Eq. jA2| )) and P (Eq. flA5] )) are seen 
to be identical. There is precisely one element with a 
(1, 2)-admissible occupation configuration in the Hilbcrt 
spaces of A and B — |2, 0) and |6, 4) respectively. Thus, 
P = (1) and the three matrices P, P and C are all of 
rank one. 

At Lf = 3, P = (—2), and is not a sub-matrix of C 
(Eq. (A3)). The matrix elements Pn and Cn are the 
coefficients of |6, 3, 3,0) and |6,3, 2, 1) in the wavefunc- 
tion \%p). The 2-body clustering constraints (14) at (3 = 3 
relate these coefficients by: 

(d 3 d + d 2 d 1 )\^) =0 

=*(Pii + Cii)|6,3) = 0. 



This relation between the single element in P and C 
proves that they have the same rank. 

The Lf = 4 case is interesting. Here P is: 



|4,0) 
13,1) 



16,2) 
1 
2 



|5,3) 
2 

-4 



P and C share the column index |5,3), but have no row 
index configurations in common. A single relation be- 
tween the row vectors of P and the row labeled by the 
partition [2, 2] in C is provided by the 2-body clustering 
conditions at /3 = 4. Without another relation, we can- 
not relate the ranks of P and C at Lf = 4. Our proof 
establishing the equality of ranks of the PEM and the 
OEM should not and is not applicable at this angular 
momentum, as Kp a = 2 and A^ = 1 with /2q = [4,0]. 



Appendix B: Rank of P 



1. (k, 2)-clustering states 



The matrices P and M were defined in Sec. llV dl as the 
particle entanglement matrices with labels (N^,Lf) in 
the (fc, 2)-admissible occupation configuration basis and 
the Jack basis. In Sec. IIV CI we showed that the PEM 



and M have the same counting; in this appendix, we 
show that P and M have the same rank. This proves 
that the counting of the PEM equals the rank of P. 

Suppose we are able to show that P = DMD', where 
D T and D' are square triangular matrices with l's on 
the diagonal and as such they have nonzero determi- 
nant. A theorem in linear algebra states that pre/post 
multiplying a matrix by one of triangular form with 
nonzero determinant leaves its rank unchanged. Thus, 
we only need to prove that P = DMD' to conclude that 
rank(P)=rank(M). 

The row and column dimensions of P and M are identi- 
cal because every (fc, 2)-admissible partition /i labels the 
Jack J". We may use partial ordering by dominance 
to order the (fc, 2)-admissible row and column configura- 
tions such that if p,^ > fa, then k < i. 

Consider a particular (k, 2)-admissible partition fii 
labeling the ith row (jth column) of P and M. Let the 
coefficient of in \J^ k ) be D^ and the coefficient of 
\uj) in \ J") be D'y. The partial ordering implies that: 



T> lk =0 if k > i 
D l2 = 1 

D'y = if I > j 

&ij = 1. 



(Bl) 
(B2) 
(B3) 
(B4) 



In other normalizations of Jack polynomials, D^ is not 
necessarily one, but is always non-zero. By the definition 
of a matrix with row-echelon form, D T and D' in row- 
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echelon form. Recall that: 

£My|J£) ® |J?> = ^Pyk) ® IK,") (B5) 

in every block of the full PEM. and |ia,) are the gen- 
eral occupation-basis states, not just the (k, 2)-admissible 
configurations. P is the sub-matrix of P labelled by 
(k, 2)-admissible partitions; therefore: 



p y =x; M «<AiJ^)^i^> (Be) 

k,l 

i.e. Py = DifcMjyD'y 



A'./ 

P = DMD' , 



(B7) 



which proves our statement that the rank of the PEM is 
given by the rank of the matrix of the coefficients indexed 
by the (As, 2)-admissible partitions. 



(k, 2)-clustering states multiplied by Jastrow 
factors 



We consider the PEM of states that are (k, 2)- 
clustering polynomials multiplied by Jastrow factors, 



rw 



Zj) M , see Eq. ([31]) . Here we identify P, a 
submatrix of the PEM with the same rank, and find that 
it contains only rows (columns) that are labeled by par- 
titions 7, (fjj) of Na (Nb) particles and angular momen- 
tum Lf {Lf) that obey the generalized Pauli principle: 
there is no more than one particle in M consecutive or- 
bitals and no more than k particles in Mk + 2 consecutive 
orbitals. The total flux N$ of the partitions ji, fjj is equal 
to the total flux of the ground state l^) being cut. 

Instead of expanding in terms of monomials, we can 
choose a different basis that incorporates all the vanishing 
properties of the Na (Nb) particles among themselves: 



k<k' 
\ k,k'<=A 

( 



) 



j% n (*i-*') m > ( Bg ) 



V 



i,i' ei 



where the Jastrow factors include only particles in part 
A and B, respectively, fa (pj) are (k, 2)-admissible 
partitions of Na (Nb) particles with angular momen- 
tum Lf (Lf ) in 2(N - 1) + MN B + 1 (for fa) and 
2(N - 1) + MNa + 1 (for fij) orbitals. The matrix 
M = (My) has the same rank as the PEM. 

Let us, for simplicity, focus on the basis states labeling 
the rows of M. The Jastrow factor can be written as 
a (1, M)-clustering Jack polynomial. Hence, both the 
Jack and the Jastrow state obey a dominance property. 



They have a root configuration with coefficient one that 
dominates any other configuration in the expansion in 
terms of occupation number states. This implies that 
also their product: 



1- 

J jXi 



no*-**')" 



(B9) 



k<k' 



has a root configuration 7; with expansion coefficient 1, 
where (7^ = (fa) 3 + M(N - j). Note that the %'s 
are precisely the configurations that label the rows of 
P. In addition, the partitions % have the same par- 
tial ordering as the //j, ie. if fa < fa, then 7^ < 7,-. 
Thus, all arguments from the previous section are ap- 
plicable here as well: There are row-echelon matrices 
D T and D' such that P = DMD', which proves that 
rank(P)=rank(M)=rank(P). 



Appendix C: The model Hamiltonian expressed as 
clustering operators 

We re- write the rotationally invariant, 2-body Hal- 
dane pseudopotential Hamiltonian, whose zero modes are 
(1, 2)-clustering states, in terms of the clustering opera- 
tors introduced in Sec. [V] Recall that the single-particle 
orbitals in the lowest Landau level form the multiplet of 
L = Na/2. In the L z basis, any 2-body interaction V 
can be expanded as: 

H = J2( m ^ N 4>/^ m 2,N 4> /2\V\m 3 ,N^/2;m 4 ,N lP /2) 

mi 

C ln 1 C m 2 Cm 3 C rn,4- (CI) 

cJ Jli is the creation operator of a single particle state of 
L z = mi,L = N&/2; 4J0) = \nn,N+/2). We first 
change basis as follows: 

\m 1 ,N <p /2- 1 m 2 ,N 4> /2) = ]T^ |mi + m 2 J) X 
x (mi +m2,£\mi,N < t > /2;m 2 ,N l i,/2) 

where (mi + m 2 , l\m\, N^/2; m%, N^/2) in the RHS are 
Clebsch-Gordon coefficients. For brevity of notation, we 
drop the label N^/2 in the superscript. To determine the 
components of the model pseudopotential in the new ba- 
sis, we recall that V is rotationally invariant (commutes 
with \L \ and L z ) and penalizes only the relative angular 
momentum of 0, thus: 



The hamiltonian in Eq. (CI I can therefore be written as: 



x (j3, N^ms, Nf/2; p - m 3 , N^/2) x 

X C mi C /3-mi C ™3 C ^-™3- 



(C2) 
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The creation and annihilation operators above create 
and destroy particles in the normalized single-particle 
orbitals. Let us denote the normalization of the single- 
particle orbital with L z — m by Af(m) . To move to the 
unnormalizcd basis, we make the transformation: 



: Af(m)c n 



(C3) 



This set of operators is identical to the 'deletion' opera- 
tors defined in Sec. VA In spinor coordinates (u,v), the 
wavefunction of the unnormalizcd orbital is: 



{u,v 



u N^/2+m v N^/2-rn^ 



(C4) 



The Clebsch-Gordan coefficients appearing in Eq. (C2) 
have the form: 

{P, N^rm, AV 2; (3 - m u = 

= KN( mi )N{(3 - miJv^VWW)!. (C5) 

K is independent of P and mi: 



K 



in 



rV4 



m^2 N *y/{N^- 1/2)! 



(C6) 



Substituting Eq. (C5) and (C3) in Eq. (C2| 



jv,. 



/3)! 



h= E E K ^ 

f)=—Nj, mx,m,3 

(A + /3)!<4_ mi d m3 Vr 



(C7) 



Comparing the equation above with the one in the text 
Eq. pL we see that f(p) = (A - /3)!(A + /?)!. 



Appendix D: Two examples of clustering constraints 

We write down the explicit relations imposed by the 
(fe + l)-body clustering constraints discussed in Sec.[v|on 
the coefficients of small wavefunctions at k = 1,2. Let us 
first consider an example at k = 1, i.e. the 1/2 Laughlin 
states. The clustering conditions are 2-body: 



E d /3-^i 

i=0 



0. 



for p = 0, 1 



T tot 

, . . . , Lt % . 



(Dl) 



Consider the N = 3 , = 6 wavefunction in the infinite 
plane geometry in which the number of orbitals is not 
restricted to + 1 = 5 as in the case of the sphere . 
The Hilbert space is spanned by 7 partitions, — 
. . . 7}. Their corresponding coefficients in \ip) are {hi, i — 
1...7}: 

\ij>) - h |6, 0, 0) + 6 2 |5, 1, 0) + 6 3 |4, 2, 0) + 6 4 |3, 3, 0) 
+fe 5 |4,l,l)+6 6 |3,2,l)+6 7 |2,2,2). 

The relations at j3 = 0, 1 respectively are: 

6ij6, 0,0) = => h = 
6a|5, 1,0) = 0^6 2 = 0. 



Thus, the clustering constraints assign zero weight to 
Ai and A2, which are not dominated by the root par- 
tition A 3 = [4,2,0] (n(A 3 ) = {10101}). The val- 
ues of P from 2 to 6 generate a set of 5 linearly 
dependent equations that fix 4 out of the 5 remain- 
ing coefficients. All the relations obtained are shown 
in Table U The solution in terms of the coefficient 
of the root partition 63 is {b\, 62, &3, &4, &5, be, bj} — 
{0,0,6 3 ,-2& 3 ,-2&3,26 3 ,-6&3}. 



TABLE I. Possible occupation configurations and the clus- 
tering conditions for the iV = 3,v = 1/2 Laughlin state at 
Ll ot — 6 on the infinite plane (no restriction to the number of 
orbitals. On the sphere the first two configurations have zero 
weight and the last two orbitals are missing as JV^ = 4 



lent of m M 


n(» 






Constraint 


61 


{2000001} 


P 


= 


61 = 


62 


{1100010} 


P 


= 1 


62 = 


bs 


{1010100} 


P 


= 2 


26 3 + b 5 = 




{1002000} 


P 


= 3 


64 + b 6 = 


br, 


{0200100} 


P 


= 4 


2b 6 + 2b- A + b 7 = 


b 6 


{0111000} 


P 


= 5 


65 + be = 


1>T 


{0030000} 


P 


= 6 


26 3 + b 4 = 



The bosonic Moore-Read state is the clustering poly- 
nomial at k = 2. The clustering conditions involve 3 
particles: 



E d p->- 

i,j=0 



jdidjlip) =0, for j9 = 0,1,. 



(D2) 



Consider the 6-particle wavefunction with L' ot = 12. 
Eq. (D2) for /3 = ensures that the weight of the par- 
titions"^, 4, 4, 0, 0, 0], [5, 4, 3, 0, 0, 0] . . . [12, 0, 0, 0, 0, 0] not 
dominated by [4, 4, 2, 2, 0, 0] is zero in the wavefunction. 
The number of such partitions whose coefficients are set 
to zero at /3 = is the number of partitions of 12 into at 
most 3 parts. Similarly, the constraints at /3 = 1 set the 
weights of the partitions [4, 4, 3, 1, 0, 0] , [5, 4, 2, 1, 0, 0] , . . . 
[11, 1, 0, 0, 0, 0] (the number of such partitions is the num- 
ber of partitions of 11 into at most 3 parts) in the wave- 
function to zero. The linear dependence of the set of 
constraints in Eq. (D2| is apparent in the fact that the 



coefficient of the partition [7, 4, 1, 0, 0, 0] is set to zero by 
a constraint at fj = and one at /3 = 1. The constraints 
at p = 11, 12 are also seen to give identical relations 
to those at P = 0, 1 for this example. The configura- 
tions [5, 4, 3, 0, 0, 0] . . . [12, 0, 0, 0, 0, 0] are only allowed in 
an infinite plane geometry. On the sphere, they would 
involve more orbitals than A^ + 1 = 5 existent ones and 
would not appear in the Hilbert space of the decom- 
position of the Moore-Read ground-state. The config- 
urations [4,4,4,0,0,0] and [4,4,3,1,0,0] appear on the 
sphere but, due to the same reason as on the infinite 
plane - that they are not squeezed from the root parti- 
tion - have zero weight. 

The 16 partitions dominated by the root partition 
[4,4,2,2,0,0] and their corresponding coefficients in ip 
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are shown in the second and first column of Table |TT] re- 
spectively. Let us discuss the 3-body clustering at (3 — 4 
in more detail: 

3(did d + 2d 3 did + d 2 d 2 d + d 2 didi)\ip) = . (D3) 

The four terms in Eq. (jD3h individually are: 



eModolV') 
d 3 did \^} 
d 2 d 2 d \^} 

d 2 didi\ip) 



= 6i 14,2,2) +6313,3,2) 
= 6 6 |4,3,1) +6 7 |4,2,2) +& 8 |3,3,2) 
= 6i|4,4,0)+6 7 |4,3,l)+6 12 |4,2,2) 
+ 6 U |3,3,2) 

= 6 2 |4,4,0) + 6 9 |4,3,1) +6io|4,2,2) 

+ 6 14 |3,3, 2) . (D4) 



The right-hand-side of each of the four terms above is a 
linear combination of different occupation configurations 
of 3 bosons with total angular momentum L' ot — /? = 8. 
Since different occupation co nfig uration states are or- 
thogonal to each other, Eq. (D3) can only be satisfied 



if the coefficient in front of every non-interacting many- 
body state is zero. Thus, we obtain four constraints on 
the coefficients from each of the four occupation config- 
urations in Eq. (|D4|): 



|4,2,2) 
|3,3,2) 
14,3,1} 
14,4,0} 



bi + 2b? + b 12 + bio 
b 3 + 2b 8 + b n + bu 
26 6 + b 7 + 6 9 = 
bi+b 2 = 0. 



(D5) 



The last relation also arises from the clustering constraint 
at P = 2. 

All the relations imposed by the clustering constraints 
at /3 = 2, . . . , 6 are shown in Table [TTJ Although not ob- 
vious, in this case as in the previous, the dimension of 
the null space of Eq. (D2) is 1. This can be analyti- 



cally proved by realizing that the Moore-Read state is 
the densest unique ground-state of a Haldane pseudopo- 
tcntial Hamiltonian which can be written in terms of the 
clustering operators 



Appendix E: Proof of step [T] in Sec. |VI B] 

We now prove the statement in Step [I] of Section [VI B| 
Kjx < Ap implies that Kp < Ap for all (k, 2)-admissible 
partitions ft that are dominated by fty. We defined fto to 
be the partition that dominates all other (k, 2)-admissible 
partitions at given Na, Lf, (Eq. ( IS ) ) : 



n(jlo) = fcO . . . fcO x(L 

2[(N A -l)/k\ 



.010...0, 



(El) 



where < x < k is fixed by the total particle number 
being Na- We are given that Ap > Kp . The case when 
Kfi = is trivial. All (fc, 2)-admissible partitions have 
distance from the cut and at least intact unit cells; 
therefore Kp < Ap for all (k, 2)-admissible partitions ft. 



TABLE II. Possible occupation configurations and the clus- 
tering conditions for the N = 6 MR state at L\ ot = 12 





n(ji) 




{20202} 


b 2 


{12102} 


63 


{20121} 


64 


{20040} 


65 


{04002} 




{12021} 


b 7 


{11211} 


f>8 


{11130} 


69 


{03111} 


but 


{02301} 


6ll 


{10320} 


612 


{10401} 


&13 


{03030} 


&14 


{02220} 


bis 


{01410} 


bid 


{00600} 





Constraint 


13 = 2: 


61 + 
63 + 


b 2 = 
ba = 


13 = 3: 


3b 3 - 
364 - 


t- 667 + b 9 = 
f- 6bg + 613 = 




662 - 


I-65 =0 


13 = 4: 


61 + 

b 3 + 


2fe 7 + 610 + bi2 = 
2b$ + bn + 614 = 




Z05 


- O7 "T OQ — U 




2b 2 - 


1- 26 7 + &9 + 6l0 = 


13 = 5: 


2b 7 - 


f-26ia +bi4 + 6i5 =0 




2b a - 


f 2b s + 613 + 614 = 




26 3 - 


|- b 6 + 67 = 




661 - 


(- 36 2 + 6b 7 + 3b 3 + 612 = 




662 - 


f 36 5 + 6b g + 3b G + 610 = 


P = 6: 


663 - 


1- 3b 6 + 66 g + 3fe 4 + 611 = 




667 - 


(- 36 9 + 6614 + 3fe 8 + 615 = 




6612 


+ 3&io + 6bi 5 + 3&u + 6i 6 = 



When Kfo > 0, we prove the required statement for 
all (k, 2)-admissible partitions by showing that Kp < Ap 
for all (fc, 2)-admissible partitions ft < fto at every Kp. 
Let us construct the partition /z (not necessarily (fc, 2)- 
admissible) at the given distance Kp = > that 
is dominated by all (fc, 2)-admissible partitions ft at the 
same distance. Assume that the orbital to the left of the 
cut is unoccupied, i.e. ni A -i(fto) = 0. If the number of 
particles to the right of the cut in fto, N r (fto), is equal 
to one then the occupation number configuration of /i is 
given by: 



n(ji) = fcO fcO (fc — 1)1 . . . (fc - 1)1 x0_^_0 I ... 01 , 

(E2) 



2 ( a mo" A m) 



where we denote the orbital cut by ' | ' in the occupation 
configuration. For N r (fto) > 1, n(fi) is: 



n(fi) = fc0.. .fc0(fc- 1)1. ..(fc - l)l(fc- 1)0 I X ... X, 

2AfI 2(A Ao -A M ) 

(E3) 



where the sequence X ... X denotes the occupation con- 
figuration of (N r (ft ) + 1) particles at distance that 
is maximally squeezed. 

As compared to n(fto), the z-angular momentum of 
the particles to the left of the cut in n(/i) is increased by 



Afj,, while that of the particles to the right of the 



cut is reduced by Kp n — Kp. Since n(p) has the same 
total z-angular momentum as n(/to): 



A Ao - A,, 
A Ao > Kp a 



K„ 



Af, > Kp. 



As every (fc, 2)-admissible partition ft with distance Kp - 
if^ that dominates [i has at least A^ intact unit cells: 



A^ > A„, Kp 



K„ 



Ap > Kp 



(E4) 
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at every distance from the cut. 

The argument for ni A _i(fl ) ^ is identical to the one 
described above. The only difference lies in the form of 
n(it): 

n(fi) = kO . . . kO (fc - 1)1 . . . [k - 1)1 0\X . . . X, (E5) 

v — x — ' v " ' 

2A " 2(A #0 -A M ) 

where the sequence X ... X is the maximally squeezed 
configuration of x + 1 particles (for N r (p, Q ) = 1) respec- 
tively k + N r (jlo) (for N r (flo) > 1) at distance K^. 



Appendix F: Proof of step [II] in Sec. |VI B| 

1. Effect of dominance on the distance from the cut 

We show that dominance, i.e. fi > // implies that 
the distance to the cut > K^, or that squeezing 
cannot increase the distance from the cut. The property 
of dominance is defined by: 



n n 



m > m' Y ^ Y v'i 



(Fl) 



i=l i=X 



for all n < N . Recall that fa > fij for i < j, where /Zj 
and /ij are the components of the partition fi. Let us 
denote the number of particles to the right of the cut for 
any partition fj, by N r (fi). The distance from the cut 
can then be rewritten as: 

K fl = ^ n m(fi>)(m- Za + 1) 

771— 2^4 

= J2fa-l A + l). (F2) 



When comparing the total distances for two partitions, 
// and jti', there are three possibilities, N r (fi) = iV r (//'), 
iV r (/i) > iV r (/i') and N r (fi) < N r ({i'). We will discuss 
them in that order: 

• N r (jj) = N r (jt'): 

M > M' Mi > X! Mi 

i=l i=l 
N r (/J,) N r (n) 

=> (Mi - *A + 1) > ^ (Mi - Ja + 1) (F3) 



Thus, > A>. 



. JV r (|i) > iV r (//): 

iV r ( M ') 

M > M' =^ X! M» > X! Mi 

i=l i=l 

=> ^ (/i, -l A + l)> (Mi - ^ + !) • ( F4 ) 



As /Zi > for all particles to the right of the cut, 

= E^Vi - «a + 1) > E^Vi - Za + 1). 
This shows that A M > K^i . 

• N r (fj,) < N r (f/y. 

fx> fx' => M» > Mi 

i=l t=l 
Nr(/0 W r (/*') 

=>> ^ ( w - i A + l) + 2 (Mi - Za + 1) 

i=l i=Af r ( M ) + l 



= K U 



> £ 04 - z A + 1) = k. 



(F5) 



The second term must be < 0, as the particles to 
the left of the cut have angular momentum /Zj < Ia- 
It is strictly negative if at least one of the \xi for 
N r (n) < i < N r (fi') is smaller that (I a — 1). 

Thus, > Kpi for every \j! that is dominated by /z. 

2. Effect of clustering conditions 

We show that the (fc + l)-body clustering conditions 
presented in the body of the paper (Eq. (14)) relate the 



partitions /i with A M > intact unit cells and distance 
K„ > from the cut to partitions // with number of 
intact unit cells given by A^ — 1 and distance from the 
cut by K^i < K^. 

Let us consider an arbitrary partition fi with A p intact 
unit cells (2A M orbitals) and distance K^: 



n(fx) ={fcO . . . fcO x ... x | x ... x ... 0} . 



(F6) 



2A„ 



-2A U 



where we placed the orbital cut after I a orbitals. In order 
to keep the discussion general, wc denote an arbitrary oc- 
cupation number configuration by the sequence x . . . x^l 
For the orbitals to the right of the cut (with angular mo- 
mentum > I a) two examples of such configurations with 



distance from the cut, K^, are: 

{ fcO. fcO x_^x | 0_^._0 10 . . 

{ fcO. fcO x_^x |Jf„0...0}. 

2A U U-2A U 



.0} 



(F7) 
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Let us now analyze the clustering condition that in- 
volve the k particles of the rightmost intact unit cell and 
the rightmost particle to the right of the cut in the par- 
tition /i (F6). We choose f3 — 2(A fl — 1) + fix and require 
the remaining Na — (k + 1) particles to occupy the same 
orbitals as in n([i). In particular, the occupation config- 
uration of the remaining particles has A M — 1 intact unit 
cells. The clustering condition relates fj, only to partitions 
that are dominated by a partition // of the forrrPB 

n{fjf) ={k0 . . . fcO(fc - 1)1 x . . . x \x . . . xO . . . 0} (F8) 



2A„ 



where x ... x is used to indicate an occupation number 
configuration where the rightmost particle to the right of 
the cut is moved to the left by one orbital. The distance 
from the cut is reduced by one: = — 1. For our 
examples in Eq. (F7|, the dominating partition is given 
by: 

n(fi') ={fc0 . . . k0(k - 1)1 x ... x | ... 10 ... 0} (F9) 



2A u -2 



l-2A u 



K u -2 



for the configuration of the first line of Eq. (F7), and: 



n{p!) ={k0_^k0(k - 1)1 x^_5 1(^-1)0... 0} 

2A„-2 Ia-2A„ 

(F10) 



for the configuration in the second line of Eq. (F7) 



Using the results from Appendix |F 1[ we conclude that 
all partitions ji 1 ^ /a involved in the clustering condition 



have A, 



A^ — 1 intact unit cells and distance from the 



cut K^i < — 1. The (k + l)-body clustering condition 
yields one constraint on the rows corresponding to all 
the involved partitions. Thus we have shown that the 
row labeled by [i can be written as a linear combination 
of the rows labeled by partitions // with K^i < and 
one less intact unit cell. 



3. Relating PEM rows to OEM rows 

Let us assemble the results of the previous appendices 
to prove the following statement: any PEM row corre- 
sponding to a partition /i with < A^ is linearly de- 
pendent on rows belonging solely to the OEM. We prove 
this statement by induction, starting with a row parti- 
tion \x with = 1 and A M > 1. Such a row partition is 
necessarily of the form: 



{fc0...fc0 x...x 1 10 ... 0} 



(Fll) 



2A U 



U-2A„ 



Using the (k, 2) clustering constraint for /? = 2fc(A M — 
1) + I a and fixing the occupation configuration of the 
remaining N — (k + 1) particles to be: 



{k0...k00Q x...x |0...0} 



(F12) 



the row partition (Fll ) can be related to row partitions 
of the OEM only. Thus, this row can be written as a 
linear combination of the rows of the OEM, independent 
on A M as long as A M > 1. 

For the induction hypothesis, let us now assume that 
all row partitions Xj at fixed K\ j = K\ > 1 with A\ j > 
K\ can be written as linear combinations of the rows of 
the OEM. Now consider a row partition fi with = 
K x + 1 and A„ > K^. 



In Appendix F 2 we have showed 



that a clustering condition involving any of the particles 
to the right of the cut and the k particles of the rightmost 
intact unit cell (to the left of the OEM cut) relates this 
partition to the partitions \j! with K^i < and A^> — 
A M — 1. This implies that the row partition fj, is a linear 
combination of the row partitions Xj . Using the induction 
hypothesis, the row partitions Xj can be written as linear 
combinations of the rows of the OEM. Therefore, the row 
partition \x can also be written as a linear combination 
of the OEM, which proves the conjecture. 



Appendix G: More clustering conditions 

We derive the clustering conditions of two particular 
states that are uniquely defined by vanishing properties 
distinct from (k,2). It should be possible to extend the 
general ideas here to other model states. For r > 3, or 
for r — 3 and k > 2, the clustering conditions derived by 
requiring that the polynomial wavefunction dies with the 
r'th power of the difference between the coordinates of 
k+ 1 particles do not uniquely define the wavefunctiorP2l. 



1. Gaffnian state 

T he b osonic Gaffnian state is a (2, 3)-clustering 
stat d 49 * 5 ^ . It vanishes as the third power between the co- 
ordinate of a cluster of two particles and that of a third 
particle approaching the cluster: 

^{2,3)(zi, Z!, z 3 , . . . , z N ) oc (zi, 3 ) 3 , for zi, 3 ->0, (Gl) 



where we define Zjj = Zi — Zj. Therefore: 



lim (zi 3 ) Q V'(2,3)(^i,^2,^3, • ■ ■ ,zn) = 0, (G2) 

21,2,Zl,3->0 



for a — 0, 1 and 2. Exactly as we did in Sec. |VJ we 
separate the coordinates z\ , zi and Z3 from the rest and 
rewrite the Gaffnian wavefunction as: 



'^(2^){Z 1 ,...,Z N ) 



2A U 



Ja-2A„ 



h,...,h Vj =i / j =i 
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where the d;. 's are the destruction operators defined in 
Section |v| Expanding Z3 3 as: 



(*1 - (^l,3))' 3 
3=0 VJ/ 



(G4) 



and inserting (G3| into Eq. (G2|, we obtain the cluster- 
ing constraints: 



E 



dp^i 2 ^i 3 di 2 di 3 \ip) = 0, V/3 > a 



(G5) 



for a — 0, 1,2. The clustering conditions at a = are 
identical to the ones we derived for the Moore- Read state 
m Sec.[vJ as a (2, 3)-clustering state also satisfies (2, 2)- 
clustering. The set of clustering conditions at each value 
of (3 are linearly dependent; in fact, for each (3 > 2, 
the number of linearly independent clustering relations 
is N r = 2. 



2. Fermionic (k, 2)-clustering states 

The fermionic counterpart of the (fc, 2)-clustered 
bosonic state is: 

rp(zi, ...,z N ) = V>(fc,2)(*i'-- -,zn) ■ - z i) ■ (G6) 

i<j 

Let us start with the simplest example, the Laughlin 
state for k = 1. From the form of the wavefunction, 
tp = rii<j( z i — z j) 3 i we see that: 

Urn zr£V(*i> = 0, for a = 0,1, 2 (G7) 

21,2-i-O 

with Zj j = Zi — Zj. Let us introduce a fermionic deletion 
operator di that destroys a fermion in angular momentum 
orbital i, analogous to the bosonic case Eq. ( 11 1. We can 
rewrite the wavefunction as: 



t/j(zi,...,z N ) = y^z l 1 1 z l 2 2 (z3 ) ...,z N \di 1 di 2 \i/ji) (G8) 
and expand z^ 2 as: 



(zi - zi,2) 1 



= E 

3=0 



l ? \ ..h-3 



(G9) 



Inserting this expression of tp into Eq. ( G6 1 and taking 



the limit z\ 2 — > 0, the only non-vanishing contribution 



is for j = a(— 0, 1, 2) — all others vanish trivially — and 
we arrive at the clustering conditions: 







E 



dhdi 2 \ip) 



(G10) 



The condition at a = is identically zero due to the 
anti-commutation relations of the fermionic operators. 
For a = 1, we find — using /? = l\ + l^. 



= 5^W^_,d,|^>, for/3>l. 



(Gil) 



When applying the above conditions, one must account 
for the anti-commutation of the fermionic deletion oper- 
ators di. Choosing a = 2 yields clustering constraints 
that are identical to those at a = 1, up to an overall 
multiplicative constant. Thus, for the fermionic model 
state at v = 1/3 we find only one clustering condition, 



Eq. (Gil I, as in the bosonic case. 



For k > 1 a very similar picture emerges. The two- 
body clustering conditions that originates from requiring: 



lim ip(zi, 
zi,2-H) 



ZN) 







(G12) 



is equivalent to Pauli exclusion statistics. In order to 
find the relevant (k + l)-particle clustering condition, we 
divide the wavefunction by a full Jastrow factor of the 
particles Zt, . . . , Zfe+i: 

0= lim f TT zt] j i,(z 1} ...,z N ). (G13) 

\*<3 ) 

Following the same steps as in the previous subsection 
we find the clustering conditions: 



fc+i 

°= e n 

h, ...,h+i 3=1 



(G14) 



with li + I2 + . . . lk+i = P- 



In principle, one can also analyze variants of Eq. ( G13 1 , 



where not a full Jastrow factor is divided out, and derive 
clustering conditions from them. However, the resulting 
conditions are identical zero due to the anti-commuting 
operators. The only non-trivial relation is the one given 
in Eq. dG14|. 



In general, when multiplying the (k, 2)-clustering 
model state with M Jastrow factors (M > 1), we find 
[M/2J 2-body clustering constraints, and (for k > 1) 
[M/2J 3-body clustering constraints, in addition to the 
original (fe+l)-body clustering constraint from the model 
state. Thus, the total number, N c , of clustering con- 
straints is iV c = 2[M/2J + 1. 
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